Parity Effect and Tunnel Magnetoresistance of Ferromagnet / Superconductor / 
Ferromagnet Single-Electron Tunneling Transistors 
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We theoretically study the tunnel magnetoresistance(TMR) of ferromagnet / superconductor / 
ferromagnet single-electron tunneling transistors with a special attention to the parity effect. It is 
shown that in the plateau region, there is no spin accumulation in the island even at finite bias 
voltage. However, the information of the injected spin is carried by the excess electron and thus the 
TMR exists. The spin relaxation rate of the excess electron can be estimated from the TMR. We 
also show that the TMR increases with decreasing the size of the superconducting island. 

PACS numbers: 72.25.-b,73.23.Hk,74.50.+r,73.23.-b 



Single-electron tunneling(SET) transistor is a key de- 
vice of "single electronics" since the transfer of a sin- 
gle electron can be controlled by gate and bias voltages 
jjl], D . The important quantity of the SET transistor is 
the electrostatic energy of excess electrons in the island 
called "charging energy" , which has a significant effect on 
charge transport, i. e. , the Coulomb blockade(CB) [Q. 
In the CB region, sequential tunneling, where tunneling 
events in each junction occur independently, is blocked 
at T = due to the increase of the charging energy. Re- 
cently much attention has been devoted to the SET tran- 
sistor with a superconducting island and normal conduct- 
ing electrodes(N/S/N) |, || §, |, @, ||. In a supercon- 
ducting island, where Cooper pairs form the condensate, 
the addition of one extra electron costs the superconduct- 
ing gap energy A. This leads to the parity effect: physical 
properties of the system depend on the parity of the elec- 
tron number in the island. In SET transistors consisting 
of normal conducting islands and electrodes (N/N/N) the 
current depends e periodically on the gate charge. How- 
ever, the tunneling current of the N/S/N SET transistor 
is expected to be 2e periodic in the gate charge below 
the crossover temperature due to the parity effect || . The 
clear signature of 2e periodicity for the N/S/N SET tran- 
sistor was observed [T 



The research on the N/S/N SET transistor has focused 
primary on the charge degrees of freedom of electrons, 
by contrast, its spin degrees of freedom have not yet re- 
ceived much attention. However, an increasing number of 
researches on spin-electronics show that the spin of elec- 
tron offers unique possibilities for finding novel mecha- 
nisms for future spin-electronic devices || [h]] . The spin- 
polarized current injected from a ferromagnetic(F) elec- 
trode into the N or S island gives rise to a nonequilibrium 
spin density in the island. In the F/N/F SET transistor, 
the tunnel magnetoresistance (TMRHs brought about by 
spin accumulation in the island dHH H |3 . Recently 
Takahashi et al. have studied the magnetoresistive effects 
of the F/S/F double tunnel junction caused by compe- 
tition between superconductivity and spin accumulation 
jL5|. They have shown that the F/S/F double tunnel 
junction is a magnetoresistive device to control supercon- 
ductivity by applying the bias voltage or current. The 
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FIG. 1: a)Schematic diagram of a F/S/F SET transistor. The 
arrows indicate the magnetizations of the left and right elec- 
trodes, b) The gate and bias voltage diagram of a symmetric 
SET transistors. Solid lines indicate the boundaries for the 
F/S/F SET transistor. Thin dashed lines are the boundaries 
for the F/N/F SET transistor (A = 0). The dot-dashed line 
indicates a boundary of the plateau region. 



suppression of the superconductivity by the spin injec- 
tion in the double tunnel junction was observed by Dong 
et al. ||). 

In this article, we theoretically study the spin- 
dependent transport of the F/S/F SET transistors at 
zero temperature with a special attention to the par- 
ity effect. We assume that the charging energy Eq is 
larg er than A and neglect the two-electron tunneling 
|17f which becomes important in the opposite situation 
(Ec < A). In such a system, we have the special re- 
gion called "plateau region" where the tunneling current 
is dominated by the transition rate from the odd-state to 
the even-state ||. We show that in the plateau region, 
spin accumulation is forbidden by the superconducting 
gap even at finite bias voltage. However, the informa- 
tion of the injected spin is carried by the excess electron 
and the TMR exists. The spin relaxation rate of the ex- 
cess electron can be estimated from the TMR. We also 
show that the TMR increases as the size of the island 
decreases. 

We consider the F/S/F SET transistor shown in Fig. 
|l| (a). For simplicity we assume that the insulating bar- 
riers of junctions 1 (left) and 2 (right) are the same; we 
subsequently set C\=Ci = C. We also assume that the 
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left and right electrodes are made of the same material 
with the spin polarization P. 

The Hamiltonian of the SET transistor is given by 

H = H L +H I + H R + H ch +H TL +H TR , (1) 

where Hl, Hj, and Hr are the Hamiltonians of the left 
electrode, the central island, and the right electrode, re- 
spectively. The Hamiltonian H c h represents the charging 
energy and the tunneling processes are described by 



H T L{R) - T kg 



(2) 



where the subscript k indicates the wave vector in the 
left (right) electrode while q represents that in the island 
and no spin flip is assumed in the tunneling process. 

The energy change due to the forward tunneling of an 
electron with spin a through the junctions 1 and 2 are 
respectively given by 



E+(n) = (l + 2n)E c -^eV g -^eV 



E^(n) = (l-2ri)E c + 



9± 



eV a 



:eV, 



(3) 
(4) 



where n is the number of excess electrons of the initial 
state in the island, Cs = 2C + C g , and the superscripts 
± implies that the number of excess electrons in the fi- 
nal state is nil. The energy changes for the back- 
ward tunneling are given by (n) — E^ (n) + eV and 
E+(n) = E+(n)+eV. 

In ordinary tunnel junctions with negligible charging 
effect, individual electrons are unstable for the forward 
tunneling across the junction, and lower their energy by 
the bias potential (— eV/2 < 0). In SET transistors, how- 
ever, there are special regions called CB region in the 
gate and bias voltage diagram, where the charge state 
of the island with n excess electrons is stable with re- 
spect to tunneling. For the symmetric SET transistor 
with a normal conducting island, the CB regions are 
the rhombuses determined by E^~(n), E^ (n) , E^ (n) , and 
E^in) > 0. The boundaries of CB regions are indicated 
by thin dashed lines in Fig. |l| (b). 

For the superconducting island, however, it is know 
that these boundaries are modified by superconductivity. 
In order to obtain the modified CB regions, we evaluate 
the tunneling rates following Fazio and Schon Since 
the island is in the superconducting state, it is conve- 
nient to rewrite the electron operators in the island in 
terms of the quasiparticle operators by using the Bogoli- 
ubov transformations. Then the transition rates are de- 
termined by using the Golden-rule arguments || . Let us 
first consider an electron tunnels from the left electrode 
into the island with n — 0, thereby changing the electron 
number from n — to n — 1 with the rate 
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where V{E) = \E\/y/E 2 - A 2 is the normalized BCS 
density of states and f(E) is the Fermi distribution func- 
tion. The tunnel conductance of the junction is defined 
as = (A-Ke 2 /h)N^N^\T\ 2 , where the tunnel matrix 
elements Tj? q and T^ q are considered as a constant T and 

denotes the density of states of the island in the 
normal conducting state (left electrode). Since we assume 
the temperature is zero, the Fermi distribution function 
can be replaced by the step function and we have 

[ (E+(0) > - A) 

r+(o) = i Je+(o) 2 -a 2 . (6) 

Therefore, the boundary of the CB region for n — given 
by E^(0) = is lifted to the solid line determined by 
£■+(0) = -A as shown in Fig. [j] (b). The other three 
boundaries also shift outward and the CB region is de- 
termined by E+(0),E-(0),E+(0), and BJ(0) > -A. 

Next we consider an electron tunnels from the island 
with ri = 1 to the right electrode, thereby changing the 
electron number from n = 1 to n — with the rate 
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x f(E - 6») [1 - /(£')] 6(E' -E + £ 2 -(l)), 



where = (Aire 2 /h)NlN?\T\ 2 and 5fi is the shift of 
the chemical potential fixed by the constraint of one ex- 
cess electron charge in the island [||. The shift of the 
chemical potential is Sfi = A and we have 



r 2 " (i) = 4 



(25 a -(l)>A) 
(1^(1)1 < A) 



£ 2 -(l) 2 -A2 



(8) 
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where d = is the average level spacing of the is- 

land. Therefore, the boundary of the CB region for 
n = l given by E% (1) = is also moved to the solid 
line determined by E^ {!) — A as shown in Fig. |l| (b). 
The other three boundaries also shift inward and the 
CB region is determined by £^(1), E± (1), E£ (1), and 
E2O-) > A. One can easily show that the CB regions 
for the other even(odd) states are spread(squeezed) like 
that for n = 0(1). 

Adjacent to the CB region, we have a so-called 
"plateau region" where the tunneling current is domi- 
nated by the tunneling rate through one junction which 
behaves like a bottleneck of the tunneling current. In 
the plateau region the tunneling current is carried by the 
even and the odd states in the following manner. An 
electron with spin a tunnels into the island from the left 
electrode. While staying in the island, the spin of the 
electron relaxes due to the spin orbit interaction, surface 
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scattering, and/or the hyperfine contact interaction [fL8| . 
After a certain time period determined by the tunnel- 
ing rate rj (T (l), the electron tunnels out of the island. 
Therefore, no spin accumulation occurs even at the fi- 
nite bias voltage. Outside the plateau region, continuous 
quasiparticle states contribute to the tunneling current 
and the spin accumulation can exist. 

We consider the tunneling current and TMR in the 
plateau region, for example, indicated by shade in Fig. |l| 
(b), which is determined by E+(0) < -A,E+(1) > A, 
and |£^~(1)| < A. In this region, the following three 
states are available: |0) three is no excess electron in the 
island, | f) there is one up-spin excess electron in the 
island, | J.) there is one down-spin excess electron in the 
island. The transition rate from |0) to |er) is given by 



r+ = r+(o) 



e 2 Ria 



Ei(0) - A 2 



and the transition rate from \a] 



r; = r 2 - (i) 



to |0) is 
d 

e 2 R 2a ' 



(9) 



(10) 



We also introduce the spin relaxation rate rj which cor- 
responds to the transition rate between | f) and | J.). 

In order to obtain the tunneling current, we construct 
the master equation for the probabilities of states po, p-\, 
and pi , which is given by 



Pa = r T p T + T l p l - | ! , - 1 ] p u 

Pa = v ips - Pa) + r+po - r-pv 



(11) 

(12) 



with the normalization condition Po + Pi + P[ = 1 • We 
calculate the stationary probabilities by requiring p a — 
Pi — Pi = 0- The solutions are 



Po 



p+pr 
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(13) 
(14) 



where 7 = E CT=n ^ ( 2I ^ + F") + T+rr + \Y~Y^ an d 
the subscript a represents the spin direction opposite to 
(7. The similar result has been obtained for the single 
discrete level system jl^]. Knowing the stationary prob- 
abilities we can calculate the tunneling current through 
the left junction 



-e Po (r+ + r+) 



(15) 



The tunneling current of the F/S/F SET transistor 
given by Eq. ( |lq ) depends strongly on whether the mag- 
netizations in ferromagnetic electrodes are parallel or an- 
tiparallel. For the ferromagnetic(F)-alignment, where 
the magnetizations are parallel, the junction resistances 
can be expressed as = i?2f = Rm and R\[ = i?2i = 
R m . Here Rm is the junction resistance for electrons in 
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FIG. 2: TMR of the F/S/F SET with (3 = 3.8 x 10" 6 is 
plotted against the normalized spin relaxation rate From 
top to bottom: P = 1.0, 0.8, 0.6, 0.4, and 0.2. 



the majority spin band and R m = (1 — P)/(l + P)Rm 
is for those in the minority spin band. We write the 
transition rates bv using the subscripts M and m as 



1 T ~ M 



and 



r r = 



p = r M /r+ ,c = 7,/r 

for the F-alignment is expressed by 
C/3{a + 1) + a/3 



T^j. Introducing a = R m /RM, 
M , the probability of the state |0) 
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C(a + l)(/3 + 2) + a(/3 + 2) (3 + 2 



(16) 



The probability Pq is independent of the spin relaxation 
rate. For the antiferromagnetic(A)-alignment, where the 
magnetizations are antiparallel, the probability of the 
even state is 



CP(a + 1) + a(3 



C(a +!)(/? + 2) + a 2 + a/3 + 1' 



(17) 



The TMR is calculated by the formula TMR. 
I a/ If- From Eqs. flf|)- (0), we have 



TMR : 



{a -If 



C(a + 1)((3 + 2) + a 2 + af3 + 1 



(18) 



The value [3 is much smaller than unity in almost all 
area of the plateau region since d -C Ec- If we assume 
that A=0.18 meV (Aluminum), £; c =8.0meV and d = 
3.0 x 10~ 5 meV [|, then we have (3 = 3.8 x 10~ 6 at the 
center of the plateau region. In Figure ||, the TMR for 
(3 = 3.8 x 10~ 6 is plotted against the normalized spin 
relaxation rate £. The TMR is monotonically decreasing 
function of £ and the spin relaxation rate of an excess 
electron rj can be estimated from the TMR by fitting the 
experimental data. If there is no spin relaxation process 
in the island, £ = 0, the TMR is approximately given by 
TMR~2P 2 /{1 + P 2 ). 

The fact that the TMR depends on the spin relaxation 
rate rj via its normalized value £ means that how much 
the spin information is transmitted is determined by the 
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FIG. 3: TMR of the F/S/F SET with P = 0.45, which is the 
typical value for Co |E0j , is plotted against the spin relaxation 
rate n. From top to bottom: the length of the island is L= 
20, 200, and 2000 nm. The area of the junction is fixed at 
20nm x 70 nmQ and the junction resistance Rm = 1 Mf2 
for L=2000 nm. We assume that A=0.18 meV (Aluminum), 
.Ec=8.0meV and d — 3.0 x 10~ 5 meV and the working point 
is set to the center of the plateau region. 



competition between the spin relaxation rate and the tun- 
neling rate through the right junction. In the plateau re- 
gion the inverse of the tunneling rate T^,* describes 
how long the excess electron with majority(minority) 
spin stays in the island. Therefore, the normalized spin 
relaxation rate £ represents the probability that the elec- 
tron with the majority spin tunnels out of the island hold- 
ing its spin direction and the TMR is a function of (. 

In the WKB approximation, the value \T\ 2 is inversely 
proportional to the length of the island L |2l], ^2|. If 
the junction parameters other than L are kept fixed, the 
density of states N^(= 1/d) is proportional to L. In this 
situation, the tunneling rate through the left junction 



does not depend on L since the size dependences of |T| 2 
and in R\ cancel out. On the contrary, the tunneling 
rate through the right junction is inversely proportional 
to L. The normalized spin relaxation rate £ decreases and 
therefore the TMR increases with decreasing L. In Fig. 
| the TMR of the F/S/F SET transistor with P = 0.45, 
which is the typical value for Co ]^(J, is plotted against 
the spin relaxation rate rj for various values of L. We 
assume that the area of the junction is fixed at 20 nm x 
70 nm and L=2000 nm, 200, and 20 nm corresponding to 
the Aluminum island with the average level spacing d = 
3.0 x 10~ 5 , 3.0 x 10~ 4 , and 3.0 x 10~ 3 meV, respectively 
Q. For the spin relaxation rate r\ = 10 7 s _1 , which is of 
the same order as that caused by the hyperfine contact 
interaction ff|, the TMR is 0.26, 2.4, and 15 % for L = 
2000, 200, and 20 nm, respectively. 

Very recently, Chen et al. reported the experimen- 
tal evidence for suppression of superconductivity by spin 
imbalance in Co/Al/Co SET transistors^). Although 
they observed the spin imbalance outside the plateau re- 
gion, their experimental results show that the effect we 
proposed is relevant to current experiments. 

In conclusion, we theoretically study the TMR of 
F/S/F SET transistors with E c > A. We show that 
in the plateau region, there is no spin accumulation even 
at the finite bias voltage. However, the information of 
the injected spin is carried by the excess electron and 
the TMR exists. The spin relaxation rate of the excess 
electron can be estimated from the TMR. We also found 
that the TMR increases with decreasing the size of the 
superconducting island. 
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